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Abstract 
Q. 

r \ , This letter suggests a new way to investigate 3-D chaos in spatial and frequency domains simul- 

Ch \ taneously. After spatially decomposing the Lorenz attractor into two separate scrolls with peaked 

spectra and a 1-D discrete-time zero-crossing series with a wide-band spectrum, it is found that 
*vq . the Lorenz chaotic attractor has an inherent frequency uniquely determined by the three system 

> ; 

parameters. This result implies that chaos in the Lorenz attractor is mainly exhibited when the 
(^ ' trajectory crosses from one scroll to another, not within the two scrolls. This is also true for 

^: 

^D ■ some other double-scroll Lorenz-like chaotic attractors, such as Chua's attractor. Some possible 

^: 

^^ \ applications of the inherent frequency and the spatial decomposition are also discussed. 

^ \ PACS numbers: 05.45.Ac, 05.45. Gg, 05.45.Pq, 05.45.Vx 
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Since the discovery of the first chaotic attractor by Lorenz [1|, 3-D chaotic attractors have 



a, 



been extensively studied to clari 
understand the essence of chaos 



how chaos occurs in the real world and to help better 
. This letter suggests a novel way to investigate chaotic 
attractors in frequency domain via a spatial decomposition method. It is found that an 
inherent frequency, i.e., a prominent spectrum peak, exists in all of the three variables of the 
Lorenz attractor (but is particularly distinguished in the 2;- variable) , and that the Lorenz at- 
tractor can be spatially decomposed into two single-scroll sub-attractors with peaked spectra 
at the inherent frequency and a 1-D zero-crossing time series with wide-band spectrum. This 
result strongly implies that, from the spectral point of view, chaos in the Lorenz attractor 
is mainly exhibited at the boundary of the two scrolls, not within them. Further numerical 
experiments show that the inherent frequency is uniquely determined by the three system 
parameters. The inherent frequency and the spatial decomposition also exist for some other 
Lorenz-like chaotic attractors with double scrolls, such as Chua's attractor y|l. Compared 

n 

with previous work on spectral analysis of "phase incoherent chaos" |^, this letter reveals 
much subtler features inherently-existing in Lorenz-like chaotic attractors. 
Consider the following Lorenz system [l|, |2| : 

X = a{y-x), (1) 

y = rx — y — xz, (2) 

z = xy — bz, (3) 

where (T,r,b are three positive parameters. When a = 10,r = 28,6 = 8/3, the Lorenz 
system has a double-scroll chaotic attractor shown in Fig. ^ Because the Lorenz system is 
invariant under the coordinate transformation [x, y, z) -^ {—x, —y, z), the chaotic attractor 
is symmetric with respect to the two planes x = and y = 0. When r > 1, the centers 
of the two scrolls, which are two equilibrium points of the Lorenz system, can be calcu- 
lated as i±^yb{r — 1), ±-\/6(r — l),r — 1 j. Roughly speaking, the double-scroll attractor 
is developed by the chaotic trajectory spirally rotating around the two points in a chaotic 
fashion. 

Although the Lorenz chaotic attractor has been extensively studied in the past three 
decades, the relationship between the chaotic behaviors and the spectra of x{t),y(t), z{t) 
has not yet been clarified ^. This letter tries to give a clearer description of this problem 
via a novel spatial decomposition method. In the following, without loss of generality, the 




FIG. 1: The Lorenz attractor 

system parameters are always set as the typical values cr = 10, r = 28, b = 8/3. 

In Fig. 121 the power spectra of x{t), y{t) and z{t) are shown (calculated by DFT with 
4-term Blackman- Harris window, the same hereinafter). It can be seen that the spectra of 
x{t) and y{t) are locally wide-band in the low frequency range and gradually approach to 
zero as the frequency increases, which shows the existence of chaos in x{t) and y(t). In 
comparison, the power spectrum of z{t) is much simpler: it has a significant spectral peak 
at the frequency fz ~ 1.3 Hz, which implies that z(t) is nearly periodic, i.e., weakly chaotic. 
Observing the time evolution of z{t) shown in Fig. El one can see that the spectrum peak 
of z{t) is a natural reflection of the nearly-invariant short-time period of z{t). 
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FIG. 2: The relative power spectra of x{t), y(t) and z{t) 



Considering the spatial symmetry of the Lorenz attractor with respect to the two planes 



time (sec) 
FIG. 3: The nearly-periodic evolution of z{t) 

X = and y = 0, the spectrum peak of z(t) implies that the rotation frequency of the chaotic 
trajectory around each scroll is nearly invariant and the mean frequency is fz- As a natural 
result, when the chaotic trajectory rotates within each scroll, the local (i.e., short-time) 
spectra of x(t) and y{t) should also have a peak at the frequency fz- Is it really true? If 
so, then it can be expected that the chaoticity of the Lorenz attractor is mainly exhibited 
when the trajectory crosses the boundary between the two scrolls. 

In the following, a simple spatial decomposition method is introduced as a tool to give 
an answer. The basic idea of the spatial decomposition method is to divide the chaotic 
trajectory into three parts: two separate single-scroll sub-attractors and a 1-D zero-crossing 
time series. The method is described as follows: 

• Two separate single-scroll sub-attractors: 

1) connect all disjoint positive parts of x{t) together to make an embedded signal 
x+(t) > 0, and connect all its disjoint negative parts to make an embedded signal 
X-{t) < (see Fig. HI for a schematic illustration on how the two embedded signals 
are extracted from the original signal x{t))] 

2) determine the sub-signals of y{t), z{t) corresponding to x+(t) and x__(t) in the time 
axis: y+(t),y^{t) and z+(t),z^(t). 

As a result, (x+(t),|/+(t), 2+(t)) constitutes a positive sub-attractor, and 
(x_(t), y_(t), z_(t)) constitutes a negative sub-attractor (see Fig. ^. Both sub- 
attractors have a single scroll. 

• A 1-D zero-crossing time series {ZC{i)}°Zi: 

1) use the zero-crossing times, i.e., the times when the chaotic trajectory crosses the 
boundary x = from one scroll to another, to create a discrete-time series {tzc{i)}'^i] 



2) ZC{i) is defined as the elapsed time between tzdi + 1) and tzcij)'- ZC{i) = tzdi + 
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FIG. 4: The extraction of a;+(i) and x^{t) from x{t) 





FIG. 5: The positive (solid) and the negative (dashed) sub-attractors of the Lorenz attr actor 



Note that the above spatial decomposition can also be made with respect to the plane 
y = to get two sub-attractors in the same way. After the above spatial decomposition of 
the Lorenz attractor, the three time series are analyzed via the spectral analysis method. 
It can be expected that the spectra of the two sub-attractors have a peak at the frequency 
of fz and the spectrum of ZC{i) is wide-band. The experimental results well support this 
prediction. Observing the power spectra of x+(t), y+(t), x_(t) and y-{t) shown in Fig. 
it can be clearly seen that a spectrum peak occurs at the same frequency fz ~ 1.3Hz in 



each spectrum. Note that the frequency fz of z{t) is actually a reflection of the nearly- 
periodic trajectory within the two sub-attractors, since z[t) is independent of the chaotic 
zero-crossing behaviors. 
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FIG. 6: The relative power spectra of x+(t), y+{t), X-.{t) and y-{t) 

Now, let us consider the power spectrum of ZC{i) shown in Fig. [7^, which keeps wide- 
band in the whole frequency range and has many spectral peaks without decaying to zero. 
Comparing all the five spectra, it is obvious that the local wide-band spectra of x{t) and 
y{t) stems from the wide-band spectrum of ZC{i). Because the elapsed time between two 
consecutive zero-crossing events is at least equal to the time of one cycle of x{t), i.e., the 
mean period of the Lorenz attractor Tz = 1/fz, the occurrence frequency of the zero-crossing 
event will not be greater than fz- It is the reason why the original spectra of x{t) and y{t) 
gradually decay to zero as the frequency increases beyond fz- 

Here, note the following fact: since the two sub-attractors are both approximately pe- 
riodic, ZC{i) will always be about n times of Tz, where n is an integer. In other words, 
ZC{i)/Tz approximates the number of rotational cycles of the chaotic trajectory within the 
current sub-attractor. This can be clearly seen from Fig. |H^. Then, one can define a new 
time series ZCn{i) = round {ZC{i)/Tz), where round(-) denotes the function rounding a real 
number to a nearest integer. In Fig. IHb, ZCn{i) is displayed to show its approximation to 
ZC{i)/Tz. Apparently, ZCn{i) has a clearer physical meaning than ZC{i), and has a similar 
spectrum to ZC{i). Figured gives a comparison of the power spectra of ZC{i) and ZC„(i), 
which are almost identical. 
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FIG. 7: The relative power spectra of ZC{i) and ZCn{i) 
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FIG. 8: ZC{i)/T, and ZCn{i) 

The above spatial decomposition of the Lorenz attractor shows that the chaotic trajec- 
tory's short-time behavior within each scroll is approximately periodic (i.e., weakly chaotic) 
and that the strongly chaotic behaviors mainly occur instantaneously when the trajectory 
crosses the boundary of the two scrolls with a mean invariant frequency not greater than f^- 
Another obvious physical meaning of the spatial decomposition is the collapse of 3-D chaos 
into 1-D space. Considering that any 3-D chaotic system has only one positive Lyapunov 
exponent, such a collapse seems apprehensible. 

A plenty of experiments were performed to verify the existence of this nearly-invariant 
spectral peak for other valid parameters of the Lorenz chaotic attractor. It is found that 
such an inherent frequency fz always exists and the spatial decomposition works well. Then, 
what is the relation between /^ and the three parameters r, a, 6? When r = 28, the surface 
of fz = F{a, b) is plotted in Fig. El As r increases, the height of the surface will rise but the 
basic shape of the surface remains. At present, an explicit mathematical formula has not 
been found to describe such a relationship. 




FIG. 9: The relation between /^ and {a, b) when r = 28 (some fluctuations may be induced by 
numerical errors) 

Is it possible to extend the above results on the Lorenz attractor to other 3-D chaotic 
attractors? For Chua's chaotic attractor generated with the following equations: 



X 

y 

z 



p{-x + y- f{x)), 

X — y + z, 

-Qy, 



(4) 
(5) 
(6) 



where f{x) = mox + 0.5(mi — mo){\x + 1| — |x — 1|), our experiments show that its inherent 
frequency also exists and the above spatial decomposition works as well. Further experiments 
have shown the existence of the inherent frequencies in some other 3-D Lorenz-like chaotic 
attractors with double scrolls. Considering prominent spectrum peaks also exist in many 
other chaotic attractors, such as the well-known phase-coherent Rossler attractors |^, y|, it 
is an interesting question to ask if the spatial decomposition-based inherent frequencies exist 
in all 3-D chaotic attractors and how to explain this delicate phenomenon theoretically. In 
other words, the question raised is: does the only one positive Lyapunov exponent mean 
"3-D chaos = 1-D chaos + 2-D near periodicity"? No definite answer is given at this time. 

Finally, some possible applications of the inherent frequency and the spatial decomposi- 
tion are discussed. 

System Identification The deterministic relationship between the inherent frequency fz 
and the system parameters reveals a new way of realizing system identification from the 
short-time waveform of any one of the three variables x{t), y{t) or z{t). This is useful for 
adaptive synchronization of chaos 
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Chaotic Cryptanalysis The identification of system parameters from the inherent fre- 
quency actually means the breaking of some chaos-based secure communication systems, if 
the system parameters serve as the secret key for encryption. In addition, another possi- 
bility of using the inherent frequency in cryptanalysis is direct extracting the plain-signal 
from the cipher-signal. For example, in chaotic modulation systems, the plain-signal is used 
to change the system parameters dynamically, which will change the short-time periods of 
the transmitted cipher-signal. By distinguishing the change of the short-time period, it is 
possible to directly estimate the plain-signal without identifying the system parameters |7[. 

Pseudo-Random Numbers Generation The wide-band spectrum of ZC{i) means that 
good pseudo-random numbers may be generated using ZC{i). A possible algorithm is to 
generate a symbolic 0-1 bit sequence according to which scroll the chaotic trajectory stays 
in J2|, Chap. 6]. Initial experiments show that the generated pseudo-random bits can pass 
some statistical tests. 

Suppressing Chaos Due to the symmetry of the two sub-attractors, the Lorenz attractor 
can merge into an approximately periodic single-scroll attractor with the following trans- 
form: x{t) -^ Gx{x(t)) and y{t) -^ Gy{y(t)), where Gx{-) and Gy{-) are even symmetric 
functions. Such a transform can fold the two scrolls from two different quadrants into the 
same quadrant and merge them into a single scroll. As a natural result, the strongly chaotic 
zero-crossing actions are suppressed and the whole attractor becomes approximately peri- 
odic at the inherent frequency fz- Apparently, the simplest transformation functions are 
Gx{x{t)) = \x{t)\ and Gy{y{t)) = \y{t)\, which correspond to the merging of the positive and 
the negative sub-attractors obtained via the spatial decomposition. 
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